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Abstract. In the effective theory of isoscalar and isovector dark matter-nucleon interactions 
mediated by a heavy spin-1 or spin-0 particle, 8 isotope-dependent nuclear response functions 
can be generated in the dark matter scattering by nuclei. We compute the 8 nuclear response 
functions for the 16 most abundant elements in the Sun, i.e. H, ^He, ^He, ^^G, 

^°Ne, ^^Na, ^^Mg, ^^Al, ^®Si, ^°Ar, ^°Ca, ®®Fe, and ^®Ni, through numerical shell model 

calculations. We use our response functions to compute the rate of dark matter capture by the 
Sun for all isoscalar and isovector dark matter-nucleon effective interactions, including several 
operators previously considered for dark matter direct detection only. We study in detail the 
dependence of the capture rate on specific dark matter-nucleon interaction operators, and on 
the different elements in the Sun. We find that a so far neglected momentum dependent dark 
matter coupling to the nuclear vector charge gives a larger contribution to the capture rate 
than the constant spin-dependent interaction commonly included in dark matter searches at 
neutrino telescopes. Our investigation lays the foundations for model independent analyses of 
dark matter induced neutrino signals from the Sun. The nuclear response functions obtained 
in this study are listed in analytic form in an appendix, ready to be used in other projects. 
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1 Introduction 

The quest for dark matter is at a turning point. Data from direct, indirect and collider 
searches for dark matter with unprecedented exposure, resolution and extension in energy 
will finally be available during the next 5-10 years [1-7]. Efficient strategies to globally 
interpret these data in terms of dark matter particle mass and interaction properties are of 
prime importance in astroparticle physics [8-10]. 

Effective theory methods have proven to be a very powerful tool in the analysis of collider 
data [11-14], dark matter direct [15-23] and indirect [24-26] detection experiments, and in 
combined studies of these different strategies [27-29]. The main advantage of the effective 
theory approach to dark matter is that it allows for a model independent interpretation of 
the different observations when all relevant interaction operators are simultaneously explored 
in multidimensional statistical analyses, as for instance in [16-18]. In contrast, comparing a 
simplistic model for dark matter to observations, important physical properties can be missed, 
and spurious correlations among physical observables can be enforced by the inappropriately 
small number of model parameters. 

In the context of effective theories for dark matter, the dark matter-nucleus interaction 
plays a special role, in that its exploration is complicated by non trivial properties related to 
the internal structure of the nuclei in analysis. The effective theory of dark matter-nucleon 
interactions [30, 31] predicts that 8 independent nuclear response functions - or form factors - 
can be generated in the dark matter scattering by nuclei. The interpretation of any dark 
matter experiment probing the dark matter-nucleus interaction is unavoidably affected by 
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the uncertainties within which the 8 nuclear response functions are known. Experiments of 
this type are dark matter direct detection experiments, and neutrino telescopes searching for 
solar neutrinos from dark matter annihilations. In this work we concentrate on the latter 
ones. 

Dark matter can be captured by the Sun while scattering in the solar medium. Dark 
matter particles accumulated in the Sun might annihilate producing a flux of potentially 
observable energetic neutrinos [32] . The solar neutrino flux from dark matter annihilations is 
strictly related to the rate of dark matter capture by the Sun (e.g. proportional to the latter, 
assuming equilibrium between capture and annihilation [32] ). It is therefore a function of the 
cross-section for dark matter-nucleus scattering, which in turn depends on the nuclear re¬ 
sponse functions computed in this work. For constant spin-independent dark matter-nucleon 
interactions, nuclear response functions for the most abundant element in the Sun are ap¬ 
proximately known [33]. For constant spin-dependent interactions, dark matter is assumed to 
scatter off Hydrogen only, and nuclei with a more complex structure are neglected. Finally, 
for momentum and velocity dependent dark matter-nucleon interactions, only simplified cal¬ 
culations have so far been performed in the literature. In Ref. [34], for instance, the rate 
of dark matter capture by the Sun is calculated for 6 momentum/velocity dependent op¬ 
erators, considering dark matter scattering from Hydrogen only. Recently, momentum and 
velocity dependent dark matter-nucleon interactions have also been explored in the context 
of helioseismology [35-37]. 

Nuclear response functions for model independent analyses of dark matter direct de¬ 
tection experiments have been calculated in [30, 38] under the assumption of one-body dark 
matter-nucleon interactions. Two-body interactions have also been included in [39, 40] in 
an investigation of spin-dependent dark matter-nucleus currents. In addition, two-body con¬ 
tributions to spin-independent dark matter-nucleon interactions have been claimed to be 
important in dark matter direct detection in [41-43]. The nuclear response functions for 
isotopes of Xe, I, Ge, Na, and F found in [30] have been applied to complementary analyses 
of current direct detection experiments [16, 20, 44], and in studies of the prospects for dark 
matter direct detection [17, 18]. 

In this paper we calculate the 8 nuclear response functions generated in the dark matter 
scattering by nuclei for the 16 most abundant elements in the Sun. We then use the novel 
response functions to calculate the rate of dark matter capture by the Sun within the general 
effective theory of isoscalar and isovector dark matter-nucleon interactions mediated by a 
heavy spin-1 or spin-0 particle. In the analysis, we comprehensively describe how the capture 
rate depends on specihc dark matter-nucleon interaction operators, and on the elements in 
the Sun. This study constitutes the first step towards robust model independent analyses of 
dark matter induced neutrino signals from the Sun. 

The paper is organized as follows. In Sec. 2 we provide the equations for computing the 
rate of dark matter capture by the Sun given an arbitrary dark matter-nucleon interaction. 
In Sec. 3 we review the effective theory of dark matter-nucleon interactions, while in Sec. 4 
we calculate the 8 nuclear response functions predicted by the theory for the most abundant 
elements in the Sun. We calculate the dark matter capture rate for all isoscalar and isovector 
dark matter-nucleon interactions in Sec. 5, and we conclude in Sec. 6. The dark matter 
response functions and the single-particle matrix elements needed in the analysis are listed 
in the Appendixes A and B, respectively. Finally, in Appendix C we provide the nuclear 
response functions of this work in analytic form. 
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2 Dark matter capture by the Sun 


Dark matter particles of the galactic halo with interactions at the electroweak scale can be 
gravitationally captured by the Sun. For a dark matter particle of mass at a distance R 
from the center of the Sun, the rate of scattering from a velocity w to a velocity less than 
the local escape velocity v{R) is given by [45] 


n: 


'{w) = TliW 0 




dai , 2 2^ 


' E)^ I Rp- 


( 2 . 1 ) 


In Eq. (2.1), Ek = m^w‘^12, dai/dE is the differential cross-section for dark matter scat¬ 
tering by nuclei of mass m* and density ni{R) in the Sun, q is the momentum transfer and 
E = (f‘/{2mi) the nuclear recoil energy. The sum in the scattering probability extends over 
the most abundant elements in the Sun, and the dimensionless parameters //j and are 
defined as follows 
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The velocity u in Eq. (2.1) is the velocity of the dark matter particle at i? —)> oo, where the 
Sun’s gravitational potential is negligible. The relation between u and wisw = + v{RY, 

and therefore VL~{w) depends on R. 

In Eq. (2.1), we consider the general case in which the differential scattering-cross section 
depends both on the momentum transfer q, and on the dark matter-nucleus relative velocity 
w. We therefore relax the assumption of constant total cross-section, commonly made in 
this context. This generalization is important in the study of arbitrary dark matter-nucleus 
interactions, as we will see in the next sections. 

Consider now a population of halo dark matter particles with speed distribution at 
infinity given by f{u). A fraction of them will be captured by the Sun, with a differential 
capture rate per unit volume given by [45] 


rdu^wn-iw). (2.3) 

dE Jo u 

The total capture rate takes the following form 

(■Rq AfJ 

C = d7Tj^ dRR^ — {R), (2.4) 

where we integrate over a sphere of radius Rq, corresponding to the volume of the Sun. The 
aim of this work is to evaluate Eq. (2.4) within the most general effective theory for dark 
matter-nucleon one-body interactions mediated by heavy spin-1 or spin-0 particles, using for 
each element in the Sun the appropriate nuclear response functions. 

From Eq. (2.4), one can calculate the differential neutrino flux from dark matter anni¬ 
hilations in the Sun. It is given by [32] 


^ = Pa. V Bf ^ 

dE^ AttD^ ^ ^ dE^ 


(2.5) 


where Ey is the neutrino energy, F^i the total dark matter annihilation rate, the branching 
ratio for dark matter pair annihilation into the final state /, and D the distance from the 
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observer to the center of the Sun. dNl/dE^ is the energy spectrum of neutrinos produced 
by dark matter annihilation into the final state /. In general, = (C'/2) tanh^(t/r), where 
t is the time variable, and r the characteristic time scale for the equilibration of dark matter 
capture and annihilation. 

In our calculations we consider the most abundant elements in the Sun, and use the 
densities ni{R) and the velocity v[R) as implemented in the darksusy code [33]. Accordingly, 
we include in the analysis the following 16 elements: H, ^He, “^He, ^‘^Ne, ^^Na, 

^^Mg, ^’^Al, ^®Si, ^^S, ^°Ar, ^°Ca, ®®Fe, and ®®Ni. Finally, we assume the so-called standard 
halo model [46], with a Maxwell-Boltzmann speed distribution for /(u), and a local standard 
of rest velocity of 220 km s“^. We leave an analysis of astrophysical uncertainties [47-49] in 
the evaluation of Eq. (2.4) for future work. 

3 Dark matter-nucleus scattering 

In this section we review the theory of dark matter scattering by nucleons and nuclei [30]. 

3.1 Dark matter-nucleon effective interactions 

The amplitude for dark matter-nucleon elastic scattering, At, is restricted by energy and 
momentum conservation, and respects Galilean invariance, i.e. the invariance under con¬ 
stant shifts of the tridimensional particle velocities. These restrictions determine how A4 
depends on the momenta of the incoming and outgoing particles. Let us denote by p (p^) 
and k (k') the initial (final) dark matter and nucleon tridimensional momenta, respectively. 
Momentum conservation implies that only three out of these four momenta are independent 
in the scattering process. A possible choice of independent momenta is p, k, and q = k — k', 
where q is the momentum transferred from the nucleon to the dark matter particle. Whereas 
q is Galilean invariant, p and k are not. Galilean invariance therefore implies that A4 must 
depend on the difference v = ^jmy^ — k/mv, rather than on p and k separately, v is the 
initial relative velocity between a dark matter particle of mass and a nucleon of mass m^. 
In addition to particle momenta, A4 can depend on the dark matter particle and nucleon 
spins, and jA, respectively. 

Any non-relativistic quantum mechanical Hamiltonian leading to a scattering amplitude 
obeying such restrictions can be expressed as a combination of the following five Hermitian 
operators 

IxAf ■ (3.1) 

The five operators in Eq. (3.1) act on the two-particle Hilbert space spanned by tensor 
products of dark matter and nucleon states, respectively |p,jx) |k, jA). The operator 
is the identity in this space, whereas and Sn denote the dark matter particle and 
nucleon spin operators. Einally, zq is the Hermitian transfer momentum operator, and v"*" 
the relative transverse velocity operator. They are Galilean invariant, and characterized by 
the matrix elements 

(p',Jx;k',jAr| *q|p,Jx;k,jiv) = (27r)^<5(k'-kp'- k-p) (3.2) 

(p',Jx;k',jjv| v-^ |p,ix;k,jiv) = (^v-k^^^e“*‘i'’'(27r)3(5(k'-bp'-k-p) (3.3) 

where is the reduced mass of the dark matter-nucleon system, and r is the position 
vector from the nucleon to the dark matter particle. Notice that energy conservation implies 
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Table 1. Non-relativistic quantum mechanical operators constructed from Eq. (3.1). Introducing the 
nucleon mass, mAr, in the equations all operators have the same mass dimension. 


V • q = —q^/{2^jq), and hence v*- • q = 0, with v*- = v -|- q/(2^Ar). This justifies the use 
of the notation v-*-. In Eqs. (3.2) and (3.3) we adopt a non-relativistic normalization for 
single-particle states. 

Only 14 linearly independent quantum mechanical operators can be constructed from 
(3.1), if we demand that they are at most linear in Stv, Sj,- and v-*-. They are listed in 
Tab. 1, and labelled as in [50], where the operator O 2 = v-*- • v-*- was neglected since it 
cannot be a leading-order operator in effective theories. They are at most quadratic in the 
momentum transfer, with the exception of O 15 , that is cubic in q. The restriction on the 
number of spin/transverse relative velocity operators is a constraint on the spin of the particle 
mediating the underlying relativistic interaction, that is assumed here to be less than or equal 
to 1. Tab. 1 also assumes that the mediating particle is heavy compared to the momentum 
transfer, i.e. long-range interactions are not included. 

The most general Hamiltonian density for dark matter-nucleon interactions mediated by 
a heavy spin-0 or spin-1 particle is hence a linear combination of 14 non-relativistic quantum 
mechanical operators. Ok, and is given by 


'H(r) 


15 r 




k=l 




(^) 



dk{r) . 


(3.4) 


In Eq. (3.4), and are the coupling constants for protons and neutrons as implemented 
in [50], and have dimension mass“^. By construction, C 2 = C 2 = 0. is the third Pauli 

matrix, and 1 denotes the 2x2 identity in isospin space. The matrices (1 ± 'r 3)/2 project 
a nucleon state into states of well defined isospin, i.e. protons and neutrons. As for the 
building blocks in Eq. (3.1), the operators Ok act on the two-particle Hilbert space spanned 
by tensor products of dark matter and nucleon states, [p, jk,jAr)) respectively. In 

the calculation of nuclear matrix elements for dark matter scattering by nuclei with well 
dehned isospin quantum numbers, it is convenient to rewrite Eq. (3.4) in terms of isoscalar 
and isovector coupling constants: 


15 

'W(r) = 

r= 0 ,l k=l 


In Eq. (3.5) d = 1, d = T 3 , and the isoscalar and isovector coupling constants, respectively, 
cl and 4 , are related to 4 = (4 + 4)/^ 4 = (4 “ 4)/2- 
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3.2 Dark matter-nucleus effective interactions 

We construct the dark matter-nucleus interaction Hamiltonian density, 7^T(r), from Eq. (3.5) 
under the assumption of one-body dark matter-nucleon interactions. Within this assumption, 
is the sum of A terms of type (3.5), one for each of the A nucleons in the target nucleus 

4 15 

^T(r) = ■ (3-6) 

i=l r=0,l k=l 


The Hermitian and Galilean invariant quantum mechanical operators 0^\r), k = 1,..., 15, 
are listed in Tab. 1. We use the index i to identify the specific nucleon to which dark matter 
couples in the scattering. Distinct nucleons are characterized by different isospin matrices 

hi)- 

Notice that within the effective theory approach reviewed here, the spin-dependent 
operators O 4 and Oq are treated independently. In contrast, the non-relativistic limit of 
a contact axial-axial dark matter-nucleon interaction generates a linear combination of the 
two operators [51]. Standard spin-dependent form factors used to interpret direct detection 
experiments account for this linear combination. The situation is different in the context of 
dark matter searches with neutrino telescopes. Here the operator is considered separately, 
in that for spin-dependent interactions dark matter is assumed to only scatter off Hydrogen 
with a constant cross-section, and obviously with no nuclear form factor. 

In the single-particle state of the ith-nucleon, it is convenient to separate the motion of 
the nucleus center of mass from the intrinsic motion (relative to the nucleus center of mass) 
of the nucleon itself. This separation induces the following coordinate space representations 
for q and v"*": 


with 


q = -f^x 5(x - y r) - i5(x -y + r)^. 


V — Vrp -|- V JY 


(3.7) 

(3.8) 



5(x — 
1 

2mv 


y + r) 
i^rS{ 




(3.9) 

(3.10) 


The operator Vx acts on the nucleus center of mass wave function at x, whereas Vy acts on 
the dark matter particle wave function at y. In Eq. (3.9), rriT is the target nucleus mass, and 
Ht the dark matter-nucleus reduced mass. Finally, the operator Vr in Eq. (3.10) acts on the 
constituent nucleon wave function at r, where r denotes the radial coordinate of the dark 
matter particle in a frame with origin at the nucleus center of mass (notice that in Sec. 3.1, r 
was the position vector from the single nucleon to the dark matter particle). Separating the 
center of mass motion from the intrinsic motion of the constituent nucleons, the only operator 
depending on the position of the nucleons relative to the nucleus centre of mass, r*, is v]^(x). 
Operators in Tab. 1 independent of v]^(x) act like the identity Ij on the ith-nncleon position 
Fj. In the coordinate space representation 1* corresponds to 5(r — r*). 


- 6 - 





Combining Eqs. (3.6), (3.8) and (3.10) with Tab. 1, we can finally write the most general 
Hamiltonian density for dark matter-nucleus interactions mediated by heavy spin-0 or spin-1 
particles as a combination of (one-body) charge and nuclear currents [30]: 


(A A 

^T(r) = X] 1 


T = 0,1 ^ 2=1 
A 


2=1 


2m AT 


• ff{i)S{r — Vi) — i 6 {r — r*) a{i) • 




2=1 
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2=1 


2m jq 


— fi) — *'5(r — r*)^r 


+ 


2 = 1 


2 mjv 


Vr X a{i) 6 {r — r^) -|- (i(r — r*) a{i) x 




(3.11) 


where a{i) denotes the set of three Pauli matrices representing the spin operator of the 
ith-nucleon in the target nucleus, and 

= c[ + i X v^y Sx cl+v^- cl + c[i 

V J rUN 


F — 
‘■M — 


loA = -2 
r_ 1 

5 “ 2 


Cy + i-A- . C14 

mN ^ 

^ ’^T 4 + Sx cl + y- -S^ cl+ Cy -b i-^ X Cg -b c[o 

m jv crijy 


-bv;i X Sv cTn -b cfq -b • Sv cC -|-^ x v;^ • Sv cE 

^ ^ mN ^ mN ^mN 

^ cl - Sx Cg 

rriN 


— Cg + Ci 2 - — X S;^ C 43 - 1—— • Sx Ci 5 
mjq ZUtv ZU/V" ZTT-tv 


(3.12) 


Inspection of Eq. (3.11) shows that dark matter couples to the constituent nucleons through 
the nuclear vector and axial charges (first line in Eq. (3.11)), the nuclear spin and con¬ 
vection currents (second line in Eq. (3.11)), and the nuclear spin-velocity current (last line 
in Eq. (3.11)). The 14 dark matter-nucleon interaction operators in Tab. 1 contribute to 
these couplings in different ways. For instance, the constant spin-independent operator Oi 
contributes to the vector charge coupling through the operator /g, whereas the constant spin- 
dependent operator O 4 contributes to the nuclear spin current coupling through the operator 

The interaction Hamiltonian relevant for nuclear matrix element calculations is finally 
obtained by integrating the Hamiltonian density (3.6) over space coordinates 


ff'j' 



(3.13) 


This latter integration eliminates the delta functions (5(r — r^) in Eq. (3.11). 


3.3 Dark matter-nucleus scattering cross-section 

From the dark matter-nucleus interaction Hamiltonian (3.13), one can calculate the ampli¬ 
tude for transitions between initial, |z), and final, |/), scattering states. We denote initial 
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nuclear states by Ikj’, J, Mj, T, Mj’), where J and T are the nuclear spin and isospin, re¬ 
spectively, and Mj and Mt the associated magnetic quantum numbers. With this notation 
|i) = Ik-r, J, Mj, T, Mt) ® |p, jxi ^x) magnetic quantum number of the dark 

matter particle, omitted in previous equations for simplicity), and an analogous expression 
applies to |/). We can therefore write 

(/l-f^Tl*) = (27r)^(5(ky p' - k-r - p) iMNR (3.14) 


with 


iMNR = 


(‘S> E 




r=0,l L i=l 

+ E ^ 

i=l 


f{i) ■ 




+ (iS) ■ E 

i=l 

(iw^E w 

2=1 
A 

(h) ■ E ^ (^: X + e-‘‘-'‘R(i) x ^„) 




2 = 1 


(3.15) 


In Eqs. (3.14) and (3.15) we use the result 

(k'r; p', jx, Mx\ i[4 v^, S;^] |kr; p, ^x) = e"*'i'''(27r)^5(k^ -b p' - kr - p) 

X (Jx>^xll[qWT,Sx] |jx>^x) > (3-16) 

and the notation ( 1 ) = (j^, M^| l[q, , S^] |j^, M^), with 1 equal to one of the operators Iq, 
I 5 , 11^, and Ig. Notice that on the right hand side of Eq. (3.16), q and = v-|-q/ {2^t) 

replace, respectively, q and in agreement with Eqs. (3.2) and (3.3). From now on |v| = w 
denotes the dark matter-nucleus relative velocity in the Sun. Importantly, each line in the 
transition amplitude (3.15) is equal to the product of a term containing information on the 
kinematics of the scattering and on the dark matter-nucleon coupling strength, i.e. ( 1 ), and 
a term given by a nuclear matrix element. 

In order to evaluate the nuclear matrix elements in Eq. (3.15), we perform a multipole 
expansion of the nuclear charges and currents using a spherical unit vector basis e\ with 
z-axis along q, and the identities 


^ V^47r(2L-b 1) i^jL{qri)YLo{^r,) 


L=o 

00 


^ A/47r(2L-b 1 ) ^'^^jL{qri)YLo{^r,) 


L=o 

00 


^ V 2 vr( 2 L + l) i 


L -2 


L =1 




XjL{qr^)YtL,{nrJ + ^x jL(gr,)Y^^,(b!rJ 


, A = ±1 
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(3.17) 


together with 

A= (A.eA)et, (3.18) 

A=0,±1 

that holds for any vector A, given a spherical unit vector basis e;^. The vector spherical har¬ 
monics in Eq. (3.17) are defined in terms of Clebsch-Gordan coefficients and scalar spherical 
harmonics: 

^LL’ii^r,) = Y.{L'mlX\L'lLM)YLUnrJex . (3.19) 

m\ 

They obey the identity = — (—1)^Y^^,^. The multipole expansion of the nuclear spin 

current, for instance, leads to 


A 

(i^) • ait) 
i=l 


Y, x/47r(2L + l)(-i)^ffi'V(g)((i^) • eo) 

L=0 

oo 

^ V27r(2L + l)(-i)^ (ASi_,;.(g) + zSi_,^,(g))((i^).e,), 

L=1 A=±l 


(3.20) 


with 




X mf^iqvi) 


i=l 
A r 
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(^)ih 




2=1 ^ 




a{i)f 


0 ) 


(3.21) 


where M^(gri) = jz,(grj)Y^^(nri) , and MiMiq^i) = jLiqri)YLMi^ri)- Assuming that 
nuclear ground states are eigenstates of P and CP, only multipoles that transform as even- 
even under P and CP contribute to the square modulus of the transition amplitude. With this 
assumption Y^M^riq) does not contribute at all, and is therefore not defined here. Expressions 
similar to Eq. (3.20) can be derived for the remaining charges and currents. Besides the two 
operators in Eq. (3.21), four additional nuclear response operators contribute to the transition 
probability, namely 


MlMiriq) 


^LMiriq) 
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i Y ( TiMiAiiq^i 
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(3.22) 
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Squaring the amplitude (3.15), summing (averaging) the result over final (initial) spin con- 
hgurations, and demanding that only multipoles transforming as even-even under P and CP 
contribute, one finally obtains [30] 


Ptotiw^,q ) = 


2jy -b 1 2 J + 1 ^ 
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(3.23) 


where the dark matter response function R^/, R^', R^', R^', R^'m^ ^ -^ae' 

are quadratic combinations of the matrix elements (1) and are defined in Appendix A. They 
depend on the momentum transfer, the dark matter-nucleus relative velocity, as well as on 
the dark matter-nucleon interaction strength. 

The nuclear response functions in Eq. (3.23) are defined as follows 


WYsiy) = AlAq) UT,Mt){J,YMt\\ BL.y{q) \\J,YMt) (3.24) 

LgSab 


where A and B correspond to pairs of operators in Eqs.(3.21) and (3.22). When A = B, 
only one letter is used. Sab = {0,2 ,...} for the pairs of operators A = B = 

A = B = r ^ ~ ^lM t’ ^ ~ ^ab = (1,3,... } for the pairs of operators 

A = B = A = B = A = B = and A = B = 

Finally, Sab = {2,4,...} for the pair of operators A = B = The integer numbers 

in Sab select multipoles transforming as even-even under P and CP. Notice that only two 
interference terms, i.e. A ^ B, can satisfy this requirement and thereby appear in Eq. (3.23). 

The nuclear response functions in Eq. (3.24) are expressed in terms of matrix elements 
reduced in the spin magnetic quantum number Mj. The reduction of a tensor operator 
Alm-,t of rank L is done by the Wigner-Eckart theorem 


{J,Mj\ALM;rAMj) = (-1)^"^-^ ^ (J| | 11 J) , (3.25) 

and it involves Wigner 3j-symbols which cancel in Eq. (3.24) after summing over spin con- 
hgurations because of their orthonormality. In the next section, we will evaluate our nuclear 
response functions using the Mathematica package of Ref. [50], which assumes the harmonic 
oscillator basis with length parameter b = y^41.467/(45A“^/^ — 25A“^/^) fm for the single¬ 
particle states. In this case, the nuclear response functions in Eq. (3.24) only depend on the 
dimensionless variable y = {bq/2Y. 
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For the rth-element in the Sun, we can finally write the dark matter-nucleus differential 
cross-section as follows 


/ 2 2 \ 


Ptot{w‘^,q^) , 


(3.26) 


which constitutes the particle physics input in the calculation of the rate of dark matter 
capture by the Sun. 


4 Nuclear matrix element calculation 

In this section we calculate the reduced nuclear matrix elements that appear in Eq. (3.24) 
for the most abundant elements in the Sun. We list analytic expressions for the associated 
nuclear response functions in Appendix C. These expressions can be used by the reader in 
analyses of dark matter induced neutrino signals from the Sun. We perform this calculation 
using the Mathematica package introduced in [50], which requires as an input the one-body 
density matrix elements (OBDME) for ground-state to ground-state transitions of the target 
nuclei in analysis. We compute these OBDME using the Nushell@MSU program [52], which 
allows for fast nuclear structure calculations based on the nuclear shell model. 

In order to relate the nuclear matrix elements in Eq. (3.24) to the underlying OBDME, 
we expand the nuclear operators in Eqs. (3.21) and (3.22), here collectively denoted by Alm;t, 
in a complete set of spherically symmetric single-particles states, |a). Here we assume the 
nuclear harmonic oscillator model for the radial part of the wave functions associated with the 
states [a). Within this assumption, single-particle states can be labelled by their principal, 
angular momentum and spin quantum numbers, respectively Uq, la and Sq, and by their 
total spin and isospin, respectively ja and ta- |a) = |uq(/qSq = lj2)jamj^-,ta = lj2,mt^). 
Here rrij^ and mt^ denote the total spin and isospin magnetic quantum numbers, whereas 
|a| represents the set of all non magnetic quantum numbers, i.e. \a) = \ \a\,mj^',mt^)■ With 
this notation, the nuclear operators in Eqs. (3.21) and (3.22) can be expanded as follows 


AlM]t 


Alm-,t \(3)al^ai3 
af} 


^ ^ (|q:| ..Axj-t-.. 

HI/31 


Kal ® ^\I3\]lM-,t 
\J{2L + 1)(2t -|- 1) 


(4.1) 


where = (-1)^^ ™i/ 3 +V 2 a|^| ^ and :: denotes a matrix element 

reduced in spin and isospin according to Eq. (3.25). The creation and annihilation operators 
Ga and ap transform as tensors under spin and isospin transformations, and their tensor 
product admits the following representation 


[af„| ® ~a\p\]LM;r = {2L + 1 ){ 2 t + 1 ) Y. 


X 


f ja L 

\—ma M mp J 


( ta r ti3 \ 
\-mt^ 0 mtj 




(4.2) 


The reduced nuclear matrix elements in Eq. (3.24) can be further reduced in nuclear isospin, 
and hence written as 


(J,T,Mr|| ALM-,r \\J,T,Mt) 


(- 1 ) 


T-Mt 


( T T T\ 

Y— Mt 0 Mx J 
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X 


(l«l 

HI/31 


{J,T:: [ 0 |q,| ® a|^|]L;T 
\/ + l)(2r + 1) 

(4.3) 


Using the definition of ground-state to ground-state OBDME, namely, 


= 


[a) I ®a\p\]L-,T 


^{ 2 L + 1 ){ 2 t + 1 ) 

we can finally write the reduced nuclear matrix elements in Eq. (3.24) as follows 


(4.4) 


{J,T,Mt\\ ALM-,r ||J,T,Mt) = (-1) 


T-Mt 


T T T 
— Mj' 0 Mx 


E 

l«ll/3| 


(|a| 


(4.5) 


which is the master equation for nuclear matrix element calculations based on the assumption 
of one-body dark matter-nucleon interactions. Since the nuclear operators Alm-,t depend on 
isospin through the matrices only, the doubly reduced matrix elements in Eq. (4.5) can 
be further simplified as follows 


{\oi\-:AL-,ry.m = \/2(2t + 1) Ml^l/2)j^\\AL\\np{lpl/2)jfi) , (4.6) 

where is the part of the operator Ax-t acting on nuclear spin and space coordinates. In 
Appendix B, we provide explicit expressions for the reduced matrix elements on the right 
hand side of Eq. (4.6), which in the case of the harmonic oscillator single-particle basis are 
known analytically, and depend on the momentum transfer through the variable y defined 
above. The Mathematica package in Ref. [50] provides an efficient implementations of these 
expressions. 

We now move on to the OBDME calculation. In this computation, the multipole number 
L is bounded from above, i.e. L < 2J, whereas r = 0,1. In contrast, the indexes ja] and 
1/31 in principle span a complete set of infinite single-particle quantum numbers. The nuclear 
shell model provides a robust framework to restrict the set of relevant jaj and |/3| in the 
OBDME definition (4.4), and to consistently truncate the infinite sums in Eq. (4.5). 

In the nuclear shell model nucleons occupy single-particle states degenerate in the total 
spin magnetic quantum number called sub-shells, or orbits. Sub-shells are solutions of the 
Schrodinger equation for a given nuclear potential (e.g. harmonic oscillator potential, Woods- 
Saxon potential, etc...) and reflect a choice of single-particle basis. Orbits are labeled with 
conventions similar to those used for atomic orbitals, e.g. the orbit Opi /2 has “principal 
quantum number” 0, orbital angular momentum 1 and total spin 1/2. Groups of energetically 
close sub-shells form major shells of progressively increasing energy. The set of fully occupied 
major shells forms the nuclear core. Eor instance, the orbits Osi/ 2 ) 0p3/2 and Opi /2 divide 
into the s and p major shells, and together form the core of, e.g, ^^Ne, which contains 8 
proton/neutron pairs. Analogously, the orbits lsi/ 2 , 0 ^ 3 / 2 , and Od ^/2 form the sd major shell, 
and the orbits lpi/ 2 , lP 3 / 2 ) O/ 5 / 2 ) and O/ 7/2 the pf major shell. Nucleons that are not in the 
nuclear core are called valence nucleons. Not all orbits are accessible to valence nucleons since 
sizable energy gaps separate adjacent major shells. Restrictions on the number of nucleons 
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Element 

2J 

2T 

P 

core-orbits 

valence-orbits 

Hamiltonian 

restrictions 

^He 

1 

1 

-b 

none 

s-p-sd-pf 

wbt [53] 

none 

^He 

0 

0 

-b 

none 

s-p-sd-pf 

wbt [53] 

none 

12C 

0 

0 

-b 

s 

P 

pewt [53] 

none 


2 

0 

-b 

s 

P 

pewt [53] 

none 

TOq 

0 

0 

-b 

none 

s-p-sd-pf 

wbt [53] 

Odg/slSi/slpO/ 

^^^Ne 

0 

0 

-b 

s-p 

sd 

w [54] 

none 

^^Na 

3 

1 

-b 

s-p 

sd 

w [54] 

none 

^4Mg 

0 

0 

-b 

s-p 

sd 

w [54] 

none 


5 

1 

-b 

s-p 

sd 

w [54] 

none 


0 

0 

-b 

s-p 

sd 

w [54] 

none 


0 

0 

-b 

s-p 

sd 

w [54] 

none 

4UAr 

0 

4 

-b 

s-p 

sd-pf 

sdpfnow [55] 

1 P 1 / 2 IP 3 / 2 O/ 5/2 


0 

0 

-b 

s-p 

sd-pf 

sdpfnow [55] 

1 P 1 / 2 IP 3 / 2 O/ 5/2 

5tiEe 

0 

4 

-b 

s-p-sd 

pf 

gxl [56] 

1 P 1 / 2 O/ 5/2 

58Ni 

0 

2 

-b 

s-p-sd 

pf 

gxl [56] 

1 P 1 / 2 O/ 5/2 


Table 2. Summary of element specific input parameters needed for the calculation of the OBDME 
via the Nushell@MSU code. We use the notation of [57] in defining the major shells. For each element 
in the Sun, we use a model space comprising the core-orbits and valence-orbits reported in this table. 
In the “restrictions” column, we list the energetic orbits not included in the calculation in order to 
make the computation numerically feasible. The interaction Hamiltonians in the next to last column 
are described in the review [57], and in the corresponding references. 


allowed in the most energetic orbits are often imposed in order to reduce the computational 
effort. The set of orbits that are actually accessible in a calculation constitutes the so-called 
model space. Therefore, the original ^-nucleon problem characterized by the bare nuclear 
interaction is simplified to a many-body problem restricted to the model space, and subject to 
an effective Hamiltonian. Effective Hamiltonians for nuclear shell model calculations can be 
computed microscopically from nuclear forces, or fitted empirically to observations. Within 
this framework, the sums in Eq. (4.5) only extend over orbits in the assumed model space, 
since the remaining states do not contribute by construction. We refer to [57, 58] for a more 
extended introduction to the nuclear shell model. 

The OBDME for orbits corresponding to the nuclear core can be analytically calculated. 
Only multipoles of nuclear response operators with L = r = 0 contribute, since in a nuclear 
core all orbits are fully occupied. One finds [59] 

'^\o^\p\ ~ \/2(2J -b 1)(2T -b l){2ja + 1) 5|q,p|(5t-o ^lo ■ (4.7) 

The calculation of the OBDME for the remaining orbits in the model space instead requires 
a numerical approach. We address this problem using the Nushell@MSU program [52, 57]. 
This code mainly relies on three inputs: the target nucleus spin, isospin and parity; the 
Hamiltonian for valence nucleon interactions (several options are provided with the code); 
the model space, including restrictions on the number of nucleons in the most energetic 
orbits. The assumptions made in our calculations are listed in Tab. 2, and closely follow 
the guidelines provided in Ref. [57] , and references therein. Assigned these inputs, the code 
first calculates the many-body ground-state wave function of the valence nucleon system 
diagonalizing the assumed interaction Hamiltonian. Then it evaluates the overlap of this 
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Figure 1. Dark matter capture rate C as a function of the dark matter particle mass for c° 7 ^ 0 
(top left panel), (top right panel), c} 7 ^ 0 (bottom left panel), and C 4 7 ^ 0 (bottom right panel). 

We report the total capture rate (thick black line), and partial capture rates specific to the 16 most 
abundant elements in the Sun. Conventions for colors and lines are those in the legends. 


wave function with the single-particle states la) according to Eq. (4.4). The OBDME that 
we obtain for 2®Si and ^®E using the Nushell@MSU w-interaction negligibly differ from 

those in the code [50] (here we use ^®E for comparison only, but it does not enter in our 
calculation). The remaining interactions in Tab. 2 were studied in [53, 55, 56]. Eor instance, 
in the full pf model space the gxl interaction was found to successfully describe binding 
energies, electro-magnetic transitions, and excitation spectra of Iron, and of various Nickel 
isotopes [56]. The major limitation of our numerical OBDME calculation hence resides in 
the use of model space restrictions. We adopt such restrictions because of limits in the 
available computing power: Nushell@MSU only runs on Windows machines, whereas our 
cluster for extensive calculations has a Unix architecture. Eor the gxl interaction, the impact 
of restrictions on observable quantitates has been discussed in [56]. For the isotopes ®®Fe and 
®®Ni, for example, a restriction of the model space where 5 or more nucleons are allowed 
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Figure 2. Left panel. Ratio of the capture rate of this work for c° ^ 0 to the capture rate computed 
with darksusy for spin-independent dark matter interactions. We report the ratio of total rates (thick 
black line), and the ratio of partial rates specific to the 16 most abundant elements in the Sun. The 
two total rates differ by at most 8 %. Right panel. Same as in the left panel, but for C 4 0. In this 
case the comparison can be performed for the total rates and for H only, since elements heavier than 
H are not included in darksusy for dark matter spin-dependent interactions. 


to be excited from the / 7/2 orbit to higher orbits implies an underestimation of the binding 
energy of the order of a few percent. 

Ultimately, the nuclear structure calculations performed here have to be considered ex¬ 
plorative, due to the restrictions in Tab. 2, and to the fact that more sophisticated interaction 
Hamiltonians could in principle be considered. At the same time, we are not aware of nuclear 
structure calculations of comparable complexity in the context of dark matter capture by the 
Sun. 

Before concluding this section, we comment on the OBDME calculation for Hydrogen. 
In the present analysis, Hydrogen constitutes a special case, in that it consists of a single 
valence nucleon system with no-core. Its OBDME can be trivially calculated as follows [59] 

'^\a\\li\ ~ ‘^I«ll7l'^ldll7l ’ 

where I 7 I corresponds to the Osi /2 orbit. 

Using the OBDME resulting from the methods outlined above, we evaluate the reduced 
nuclear matrix elements in Eq. (3.24), and hence the dark matter-nucleus scattering cross- 
section (3.26) for the most abundant elements in the Sun. This cross-section will allow us 
to calculate the rate of dark matter capture by the Sun (2.4) for all interaction operators in 
Tab. 1, as we will see next. The nuclear response functions that we obtain in this analysis, 
i.e. Eq. (3.24), are listed in Appendix C, and can be used by the reader for other projects. 

5 Numerical evaluation of the capture rate 

In this section we numerically evaluate the dark matter capture rate by the Sun, Eq. (2.4), 
using the nuclear response functions derived in the previous section, and collected in analytic 
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Figure 3. Same as in Fig. 1, 


but for the interaction operators O 3 and O 5 . 


form in Appendix C. We study one operator at the time, and for each interaction operator in 
Tab. 1, we separately consider the corresponding isoscalar and isovector coupling constants. 
In the figures, we report the dark matter capture rate as a function of the dark matter 
particle mass, varying in the range 10 - 1000 GeV. When a coupling constant is different 
from zero, it takes the reference value of 10“^ with = 246.2 GeV. Using the same 
interaction strength in all panels allows for a straightforward comparison between capture 
rates associated with different operators. For definiteness, we assume a spin = 1/2 for the 
dark matter particle. 

5.1 Constant spin-independent and spin-dependent interactions 

We start with the capture rate for the interaction operators Oi and O 4 , corresponding to 
constant, i.e. velocity and momentum independent, dark matter-nucleon interactions. Fig. 1 
shows the capture rate C as a function of for the two operators. The top panels refer to 
the couplings constants c? and C 4 , whereas the bottom panels correspond to c\ and c\. In the 
plots we report the total capture rate (thick black line), and partial capture rates specific to 
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Figure 4. Same as in Fig. 1, but for the interaction operators Og and Oj. 


the 16 most abundant elements in the Sun. Conventions for colors and lines are those in the 
legends. 

In the case 7 ^ 0 many elements contribute to C in a comparable manner. The leading 
contributions come from ^He, and ^®Fe, with an additional sizable contribution due to 
^^Ne for > 400 GeV. For C 4 7 ^ 0 the most effective isotopes in capturing dark matter are 
H and though the latter significantly contributes for rriy^ > 100 GeV only. Similarly, in 
the case c} 7 ^ 0 the most important element is H, though also ^®Fe gives a sizable contribution 
to C for large values of m^. Finally, for C 4 7 ^ 0 only H is relevant in the dark matter capture 
by the Sun. 

Fig. 2 compares the isoscalar rates of Fig. 1 with the spin-independent and spin- 
dependent capture rates computed by darksusy. For constant spin-independent interactions, 
corresponding to the Oi operator, darksusy uses a simplified version of Eq. (2.1), namely 


(w) = riiW 0 
i 




1 


dEai —^exp(- 2 ?/), 

J/w'^ 


(5.1) 
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Figure 5. Same as in Fig. 1, 


but for the interaction operators Os and Og. 


where Uj is the total dark matter-nucleus scattering cross-section in the limit of zero momen¬ 
tum transfer, y = {bq/2)‘^, and 


b 



0.91 


/ mi y/3 
VGeV/ 


-bO.3 


fm, 


(5.2) 


which allows to analytically compute Vl~{w). In the case of constant spin-dependent inter¬ 
actions, corresponding to the O 4 operator, darksusy calculates the capture rate for H only, 
and neglects other elements. Other interaction operators are not included in the program, 
and cannot be used for comparison. 

The left panel of Fig. 2 shows the ratio of the capture rate of this work for c)* 7 ^ 0 to the 
capture rate computed with darksusy for spin-independent interactions. We report the ratio 
of total capture rates, and the ratio of partial rates specific to different elements in the Sun. 
Whereas for elements like ®®Fe and ®®Ni the two rates differ up to 25% for ~ 1 TeV, the 
total rate computed with our nuclear response functions and the one obtained from Eq. (5.1) 
differ by at most 8 %. We conclude that for constant spin-independent dark matter-nucleon 
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Figure 6. Same as in Fig. 1, but for the interaction operators Oio and On. 

interactions, the capture rate is only moderately affected by the use of refined nuclear response 
functions. 

The capture rate for constant spin-dependent dark matter interactions computed with 
darksusy is systematically smaller than the capture rate of this work for ^ 0. This effect 
is however important for dark matter masses larger than 100 GeV only. Neglecting elements 
heavier than H, and in particular induces an error on the total capture rate of about 
25% for ~ 1 TeV, as shown in the right panel of Fig. 2. 

In summary, the capture rate for the operators Oi and O 4 found with the nuclear 
response functions of this work does not dramatically differ from that of previous studies. 
However, in the future errors at the 10-20% level on the capture rate induced by simplistic 
form factors might non negligibly alter the interpretation of a hypothetical signal in terms 
of dark matter particle mass and interaction properties. 
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Figure 7. Same as in Fig. 1, but for the interaction operators O 12 and O 13 . 


5.2 Velocity and momentum dependent interactions 

We now move on to our results for the capture rate of the operators Oi, i = 3, 5 ..., 15. We 
report these results in Figs. 3, 4, 5, 6, 7, and 8, which show total and partial capture rates 
as a function of the dark matter particle mass. In each panel the thick black line repre¬ 
sents the total capture rate, whereas partial rates correspond to colored lines, as explained 
in the legends. Inspection of these figures shows that the most important element in the 
determination of C significantly depends on the dark matter-nucleon interaction operator, 
on whether the coupling is of isoscalar or isovector type, and on the value of rriy^. Elements 
that contribute the most to the capture rate for at least one interaction operator, and in a 
specihc dark matter particle mass range are H, ^He, ^^Al and ®®Fe. The existence 

of a variegated sample of elements important in the dark matter capture process shows the 
significance of detailed nuclear structure calculations. This conclusion is in particular true for 
interaction operators that favor dark matter couplings to nuclei heavier, and more complex 
than H or '^He. 

The properties of the 6 nuclear response operators in Eqs. (3.21) and (3.22), and the solar 
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Figure 8 . Same as in Fig. 1, but for the interaction operators O 14 and O 15 . 


nuclear abundances in Tab. 3 determine the most important isotopes for a given operator. In 
the small momentum transfer limit the response operator measures the mass number 

A of the nucleus, and is therefore larger for heavy elements, like for instance ®®Fe. Operators 
coupling via Mlm;t are Oi, O 5 , Og, and On, though with different velocity and momentum 
suppressions. For these operators a compromise between nuclear abundance and mass number 
determines the most relevant elements in the capture process. The response operators 
and measure the nucleon spin content of the nucleus, and favor nuclei with unpaired 

protons or neutrons, like H and These isotopes are important for operators like O 4 , Oq, 
Oy , Og, Oio, Oi 3 , and £>14, that couple via Similar interpretations exist for 

the remaining nuclear response operators. For instance, Alm;t measures the nucleon angular 
momentum content of the nucleus, and the content of nucleon spin-orbit coupling in 

the nucleus [21]. It can be shown that favors nuclei with an unpaired nucleon in a non 

s-shell orbit, whereas r fetors heavy elements with orbits of large angular momentum 
not fully occupied [21]. For these reasons the element ®®Fe is particularly important for 
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Element 

Average mass fraction 

Element 

Average mass fraction 

H 

0.684 

^4Mg 

7.30x10-4 

^He 

0.298 

"Ul 

6.38x10-5 

^He 

3.75x10-4 

^«Si 

7.95x10-4 

12 c 

2.53x10“^ 

32S 

5.48x10-4 


1.56x10“^ 

4UAr 

8.04x10-5 


8.50x10“^ 


7.33x10-5 


1.92x10“^ 

5B|r 

1.42x10-5 


3.94x10“^ 


8.40x10-5 


Table 3. List of average mass fractions for the 16 most abundant elements in the Sun as implemented 
in the darksusy program [33]. The underlying solar model is introduced in [60]. 


interaction operators that generate the nuclear response operator like O 3 , O 12 , and 

Oi5- 

For elements up to we assume the solar abundances reported in Ref. [60]. For 
heavier elements we consider the abundances of Ref. [61]. These are the abundances im¬ 
plemented in the darksusy program, which we use to calculate the average mass fractions 
in Tab. 3. Capture rates linearly depend on the radial number densities rii (see Eq. (2.1)), 
which are in turn proportional to the corresponding mass fractions. Assuming a different 
solar model, i.e. different mass fractions, would impact our results accordingly. Conservative 
relative uncertainties on the solar abundances are listed in Tab. 4 of [62], and range from 
11.8% for ®®Fe to 45.3% for ^^Ne. Smaller uncertainties are quoted in [63]. Elements not 
included in Tab. 3 (heavier or lighter than ^®Ni) have abundances at least a factor of a few 
smaller than ®®Ni, and are neglected in all present calculations. Whether the corresponding 
nuclear response functions can compensate for the small abundances of these isotopes, is an 
interesting question that we leave for future work. 

Also the behavior of the capture rate as a function of the dark matter particle mass 
strongly depends on the nature of the dark matter-nucleon interaction. In the log-log planes 
of Eigs. 3, 4, 5, 6 , 7, and 8 , we observe steeply decreasing lines, e.g. for c) / 0, roughly 
flat lines, e.g. for 7 ^ 0 , bumps, e.g. for C 3 7 ^ 0 , and even more complex behaviors, like 
in the case of Cg ^ 0. Different factors intervene in determining the exact dependence of 
the capture rate on including what element dominates the capture process, its nuclear 
structure and the resulting nuclear response functions, and the intrinsic momentum/relative 
velocity dependence of the operator in analysis. 

Another important result of this work is to observe that the operators Oi and O 4 
do not necessarily dominate the dark matter capture process. We find that the operator 
Oil = -q/mjv generates a total dark matter capture rate larger than that associated with 
the operator O 4 for values of the dark matter particle mass larger than approximately 30 GeV. 
This result is clearly illustrated in Eig. 9, where we compare the total dark matter capture 
rate as a function of rriy^ for the operators Oi, O 4 and On assuming isoscalar interactions. As 
in the previous figures, we consider the same value of the coupling constant, i.e. 10 “^ 
for the three operators. The relative strength of the three interactions in Eig. 9 is hence 
determined by the matrix elements of the nuclear response operators 

and when evaluated for the most abundant elements in the Sun, and by the intrinsic 

momentum/relative velocity dependence of the three operators. Notice that the response 
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Figure 9. Total capture rate for the interaction operators Oi, O 4 , and On. In the three cases we 
assume the same value for the isoscalar coupling constant, i.e. 10“^ m“^, with rriy = 246.2 GeV (we 
set to zero the isovector coupling constant). The operator On, though never included in experimental 
analyses, generates a capture rate larger than that associated with the operator O 4 for > 30 GeV. 


operator Mlm;t{q) affects the cross-sections generated by Oi and On, whereas a linear 
combination of determines the cross-section associated with O 4 . 

6 Conclusions 

We have calculated the 8 nuclear response functions generated in the dark matter scattering 
by nuclei, i.e. Eq. (3.24), for the 16 most abundant elements in the Sun. We have carried 
out this calculation within the general effective theory of isoscalar and isovector dark matter- 
nucleon interactions mediated by a heavy spin-0 or spin-1 particle. This theory predicts 14 
isoscalar and 14 isovector dark matter-nucleon interaction operators with a non trivial de¬ 
pendence on velocity and momentum transfer. In contrast, current experimental searches for 
dark matter focus on 2 constant spin-independent and spin-dependent interaction operators 
only. 

We have used the nuclear response functions found in this work to calculate the rate 
of dark matter capture by the Sun for the 14 isoscalar and the 14 isovector dark matter- 
nucleon interactions separately. We find that different elements contribute to the dark matter 
capture rate in a significant manner. H, ^He, ^'^Al and ^®Fe generate the leading 

contribution for at least one interaction operator, and in a specific dark matter particle mass 
range. Nuclear structure calculations, like those performed in this work, are hence crucial to 
accurately compute the rate of dark matter capture by the Sun, in particular for interaction 
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operators that favor dark matter couplings to nuclei heavier, and more complex than H or 
^He. 

Another important result found in this work concerns the operator On = • q/mAr, 

which couples to the nuclear vector charge operator. For > 30 GeV, this operator 
generates a capture rate larger than the rate induced by the operator O 4 = • Sw, i-e. the 

constant spin-dependent operator commonly considered in dark matter searches at neutrino 
telescopes. This result was not known previously, and should be kept in mind in the analysis 
of dark matter induced neutrino signals from the Sun. It is however not unexpected, in that 
Oil is independent of the nucleon spin, i.e. C oc and of the transverse relative velocity 
operator. 

Our findings significantly extends previous investigations, where the dark matter capture 
rate was calculated for constant dark matter-nucleon interactions only (see however [34] for 
an interesting exception), and using simplistic nuclear form factors. The nuclear response 
functions obtained in this work are listed in analytic form in Appendix C, and can be used 
in model independent analyses of dark matter induced neutrino signals from the Sun. 
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A Dark matter response functions 

Below, we list the dark matter response functions that appear in Eq. (3.26). The notation is 
the same used in the body of the paper. 
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B Single-particle matrix elements of nuclear response operators 

Here we list the single-particle matrix elements of the nuclear response operators in Eqs. (3.21) 
and (3.22). Eqs. (B.2), (B.3), (B.4), and (B.5) are implemented in the Mathematica package 
of Ref. [50]. 

Only 4 independent nuclear response operators are actually generated in the dark 
matter-nucleus scattering, and need to be considered in order to evaluate the dark matter- 
nucleus scattering cross-section. These are MjM{qT^i), ' ^(0) 

lVIj^_i_]^(grj) • f(T(i) X This result follows from the identities 



kjM-k^) -y 



2J+ i ^36-i(grO (B.l) 


2=1 


The reduced single-particle matrix elements of the four independent nuclear response op¬ 
erators are given in the following, where to simplify the equations we use the notation 
(ajj • 11/3) = {nklkk)jo\\ ■ ||n/ 3 (// 3 l/ 2 )j^), and [A] = \/2A -f 1, for any index A. They 
read as follows 
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Eqs. (B.2), (B.3), and (B.4) also appear in the calculation of nuclear matrix elements for 
electroweak lepton-nucleus interactions. The latter expression is instead needed to evaluate 
the matrix elements of the nuclear response operators and , specific to the dark matter- 
nucleus scattering. Different combinations of Wigner 3j, 6j and 9j symbols appear in the 
equations above, which also depend on residual radial matrix elements of spherical Bessel 
functions and of their derivatives at p = grj. In the case of the harmonic oscillator single¬ 
particle basis, these radial matrix elements can be analytically evaluated as follows 
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+(2Z/3 + 2/c + l)ii^i[-(L — Ip — la — 2k — 2k + 1); L + -; y]!- , 
where y = {qb/2)‘^, and i-Fi is the confluent hypergeometric function. 


C Nuclear response functions 

Below, we only list nuclear response functions different from zero. 


Hydrogen (H) 

W^{y) = 0.0397887 
Wl}{y) = 0.0397887 
Wl^{y) = 0.0397887 
W^{y) = 0.0397887 

Helium (^He) 

W^{y) = 0.358099e-2j/ 
Wl}{y) = 0.0397887e“2^ 
Wl^{y) = 0.119366e"22' 
W^{y) = 0.119366e"22' 

Helium (^He) 

W^{y) = 0.31831e“2y 


WOO(y) = 0.0397887 
W^l{y) = 0.0397887 
W^O(y) = 0.0397887 
W^l{y) = 0.0397887 


ppOO(y) = 0.03978876-2?^ 
Wl},{y) = 0.03978876-2?? 
7^10 (y) = -0.03978876-2?? 
W^},{y) = -0.03978876-2?? 


W^{y) = 0.0795775 
Wl}{y) = 0.0795775 
W^{y) = 0.0795775 
W^{y) = 0.0795775 

(C.l) 


W^{y) = 0.07957756-2?? 
W^{y) = 0.07957756-2?? 
W^{y) = -0.07957756-2?? 
W^}{y) = -0.07957756-2?? 

(C.2) 
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Carbon (^^C) 

= 0.565882e“2^(2.25 - yf 
Wl^,{y) = 0.0480805e-2y 
W^^„{y) = e“22'(-0.371134 + 0.164948y) 

Nitrogen (^^N) 

b/^0(y) = e"22'(11.6979 - 11.1409y + 2.67574y2) 

W^,{y) = 0.0230079e"2y(l.20986 + yf 
W^{y) = 0.134532e“2?^(0.707578 - yf 
Wl^fy) = 0.0905048e“2y 
Ty|0(y) = 0.00126432e“2*' 

W'^^{y) = 0.0424075e“2y 
^MV(y) = e-22'(-1.02414 + 0.483267?/) 

W^^{y) = e“22'(0.0534451 - 0.0755325?/) 

Oxygen (^®0) 

W'^iy) = 0.000032628e“2^(395.084 - 200.042?/ + y‘^f 
Wl^fy) = 0.000032628e“2^(3.66055 - yf 
W^^fy) = e"22'(-0.0471874 + 0.0367831?/ - 0.00664641?/2 + 0.000032628?/^) 

Neon (^°Ne) 

W'^iy) = 0.0431723e“2^(13.5766 - 9.05108?/ + y'^f 
Wl^fy) = 0.003480776-22/(2.50001 - yf 
W^^fy) = e-222(-0.416077 + 0.443815?/ - 0.1416?/2 + 0.0122586?/^) 

Magnesinm (^^Mg) 

W'^iy) = 0.123467e-2?2(9.63385 - 7.49299?/ + y^f 
Wl^fy) = 0.02608166-22/(2.5 - yf 
W^^fy) = 6-222(-1.36673 + 1.6097?/ - 0.567072?/2 + 0.056747?/^) 


(C.5) 


(C.6) 


(C.7) 


(C.8) 
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Sodium (^^Na) 

Ty“(y) = e"22'(42.0965 - 63.4498y + 32.5913^^ - 6.57878/ + 0.483166/) 

W}^{y) = e"22'(0.0795776 - 0.212207?/ + 0.182941/ - 0.0543892/ + 0.00523012/) 
W}^{y) = e“22'(-1.83028 + 3.81972y - 2.50445/ + 0.597822/ - 0.04545/) 

W'^iy) = e“22'(-1.83028 + 3.81972y - 2.50445/ + 0.597822/ - 0.04545/) 
irOO(y) = e“22'(0.0126672 - 0.0262533?/ + 0.0401886/ - 0.010514/ + 0.00078605/) 
W^,{y) = e“2*'(0.00917577 - 0.0167053?/ + 0.0332751/ - 0.00765719/ + 0.000597676/) 
W^,{y) = e"22'(0.0107811 - 0.020986?/ + 0.0360971/ - 0.00876213/ + 0.000626718/) 
W^l{y) = e"22'(0.0107811 - 0.020986?/ + 0.0360971/ - 0.00876213/ + 0.000626718/) 
W^{y) = e"22'(0.0253345 - 0.0750847?/ + 0.100235/ - 0.0384261/ + 0.00466396/) 

W^{y) = e“22'(0.0183515 - 0.0567009?/ + 0.0887794/ - 0.0374699/ + 0.00477955/) 

W^{y) = e“2*'(0.0215622 - 0.0652627?/ + 0.0941439/ - 0.0379511/ + 0.00472138/) 

W^}{y) = e“22'(0.0215622 - 0.0652627?/ + 0.0941439/ - 0.0379511/ + 0.00472138/) 

Wl^,{y) = e"22'(0.612149 - 0.49308?/ + 0.107832/) 

Wll{y) = e"22'(0.00940911 - 0.00747826?/ + 0.00163204?/2) 

Wl^,{y) = e"22'(-0.075893 + 0.060682?/ - 0.0110124/) 

Wll{y) = e"22'(-0.075893 + 0.060682?/ - 0.0110124/) 

Wf,{y) = e“22'(0.000495589 - 0.00010394?/ + 0.00000544981?/2) 

Wf,{y) = e“22'(0.00000616583 + 0.00008381?/ + 0.0002848/) 

Wf,{y) = e“2*'(-0.0000552785 - 0.000369894?/ + 0.0000393968?/2) 

Wf,{y) = e“2*'(-0.0000552785 - 0.000369894?/ + 0.0000393968?/2) 
lyO/y) = e“22'(0.0335711 - 0.0268568?/ + 0.00656896/) 

W^{y) = e“22'(0.00772326 - 0.00617861?/ + 0.0021619/) 

W^{y) = e“22'(0.0161021 - 0.0128817?/ + 0.00362952/) 

Wf{y) = e"22'(0.0161021 - 0.0128817?/ + 0.00362952/) 

W^^„{y) = e"22'(-5.07498 + 5.86765?/ - 2.09908/ + 0.226345/) 

= e"22'(-0.0273574 + 0.0474719?/ - 0.0213121/ + 0.00280825/) 

W]^^„{y) = e“22'(0.220651 - 0.382932?/ + 0.17682/ - 0.0226015/) 

W^^„{y) = e“2*'(0.62922 - 0.727336?/ + 0.243236?/2 - 0.0210943/) 

W^^{y) = e“22'(-0.0291634 + 0.0548817?/ - 0.0305345/ + 0.00476387/) 

^#a(/ = e"22'(-0.0119052 + 0.0231539?/ - 0.0164035/ + 0.00310235/) 

W^^{y) = e"22'(-0.024821 + 0.0482732?/ - 0.02884?/^ + 0.00481368/) 
lyoi^(y) = e"22'(-0.013988 + 0.0263236?/ - 0.0171362/ + 0.00306717/) 

(C.9) 
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Aluminium (^^Al) 


wOO 
w M 

(y) 


(y) 
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yy M 

(y) 

wOl 

w M 

(y) 
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U/OO 
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(y) 

wlO 
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= e"22'(87.0146 - 146.097y + 83.5367^^ - 18.5981y^ + 1.43446/) 

= e"22'(0.119366 - 0.31831y + 0.337291/ - 0.132526/ + 0.018155/) 

= e“22'(-3.22283 + 7.00266y - 4.92756/ + 1.33587/ - 0.11524/) 

= e“22'(-3.22283 + 7.00266y - 4.92756/ + 1.33587/ - 0.11524/) 

= e“22'(0.0309465 - 0.0367242y + 0.0265347/ - 0.00241606/ + O.OllOOll/) 

= e“2*'(0.0218834 - 0.00944476y + 0.011506/ + 0.000953537/ + 0.0104813/) 

= e"22'(0.0260233 - 0.0210567y + 0.0158643/ + 0.000606077/ + 0.0105713/) 

= e"22'(0.0260233 - 0.0210567y + 0.0158643/ + 0.000606077/ + 0.0105713/) 

= e"22'(0.0618929 - 0.210848?/ + 0.244466/ - 0.0942682/ + 0.0243737/) 

= e“22'(0.0437667 - 0.165622?/ + 0.221193/ - 0.101991/ + 0.0277477/) 

= e“22'(0.0520466 - 0.18713y + 0.233007/ - 0.0985082/ + 0.0259327/) 

= e“22'(0.0520466 - 0.18713y + 0.233007/ - 0.0985082/ + 0.0259327/) 

= e"22'(2.80498 - 2.24306?/ + 0.455491/) 

= e"22'(0.021493 - 0.0156159?/ + 0.00596886/) 

= e"22'(-0.180417 + 0.137389?/ - 0.0239615/) 

= e"22'(-0.180417 + 0.137389?/ - 0.0239615/) 

= e“22'(0.0000680703 - 0.000376682?/ + 0.00340251/) 

= e“22'(0.0149622 - 0.00563307y + 0.00440385/) 

= e“2*'(-0.0010092 + 0.00298228?/ + 0.00281525/) 

= e“2*'(-0.0010092 + 0.00298228?/ + 0.00281525/) 

= e“22'(0.126043 - 0.100835y + 0.0237577/) 

= e“2*'(0.05736 - 0.045888y + 0.012102/) 

= e“22'(0.0850285 - 0.0680228?/ + 0.016845/) 

= e"22'(0.0850285 - 0.0680228?/ + 0.016845?/2) 

= e"22'(-15.6228 + 19.3589?/ - 7.23234/ + 0.79705/) 

= e"22'(-0.0370794 + 0.0852545?/ - 0.0449284/ + 0.00866992/) 

= e“22'(0.578632 - 1.00438y + 0.491252?/2 - 0.0730693/) 

= e“22'(1.00112 - 1.15934?/ + 0.40275/ - 0.0364952/) 

= e“22'(-0.0883243 + 0.185775y - 0.104001/ + 0.0163635/) 

= e"22'(-0.0501045 + 0.114845y - 0.0729898/ + 0.0131315/) 

= e"22'(-0.0742731 + 0.170242y - 0.105744/ + 0.0188197/) 

= e"22'(-0.0595834 + 0.125323y - 0.0717204/ + 0.011398/) 

(C.IO) 
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Silicon (^®Si) 


W^{y) = 0.281695e"2^(7.44089 - 6.37784y + 7/2)2 
Wl^,{y) = 0.0739103e-2y(2.5 - yf 
W^^„{y) = e“22/(-2.68415 + 3.37434y - 1.281y2 + 0.1442927/^) 

Sulfur (32S) 

Ty^O(y) = 0.580305e“22^(5.92494 - 5.43118y + y'^f 
Wl^,{y) = 0.0765941e-2y(2.5 - yf 
W^^„{y) = e“22/(-3.l2284 + 4.11173y - 1.67217/2 + 0.2108277/^) 

Argon (^°Ar) 


(C.ll) 


(C.12) 


W^{y) = e"2s/(3l.8294 - 65.96187/ + 48.5834y2 _ 15.1947/^ + 1.9036/ - 0.0595886//^ 

+ 0.000544329/) 

Wl^{y) = e“22/(0.318304 - 1.06524y + 1.24846/ - 0.62249/ + 0.141618/ - 0.0138797/ 
+ 0.000480513/) 

W}f{y) = e“22/(-3.l8299 + 8.62425y - 8.02539/ + 3.19316/ - 0.554467/ + 0.0353797/ 
-0.000511426/) 

Wf^{y) = e"22/(-3.l8299 + 8.62425y - 8.02539/ + 3.19316/ - 0.554467/ + 0.0353797/ 

- 0.000511426/) 

Wl^fy) = e"22/(0.299629 - 0.373798// + 0.154895/ - 0.0238983/ + 0.00122474/) 

Wll{y) = e“22/(0.00414999 - 0.0181474y + 0.0240755/ - 0.00926264/ + 0.00108115/) 
Wl^fy) = e“22/(-0.0352627 + 0.0990955// - 0.0683453/ + 0.0161561/ - 0.00115071/) 
Wll{y) = e“2s/(-0.0352627 + 0.0990955// - 0.0683453/ + 0.0161561/ - 0.00115071/) 
Wf}^„{y) = e“2j/(-3.08821 + 5.12625// - 2.89248/ + 0.653386/ - 0.0526576/ 

+ 0.000816493/) 

W}}^„{y) = e"22/(-0.036345 + 0.140282// - 0.1719177/2 + 0.0770456/ - 0.0134973/ 

+ 0.000720769/) 

W]f^„{y) = e“22/(0.308826 - 0.709394// + 0.515378/ - 0.153134/ + 0.0185641/ 

- 0.000767139/) 

W^^fy) = e“22/(0.363444 - 1.17124// + 1.09117/ - 0.373592/ + 0.0452762/ 

- 0.000767139/) 

(C.13) 
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Calcium (^°Ca) 


W^{y) = 0.000016743e“2^(1378.8 - 1387.5% + 281.953^^ - 
irOO(y) = 0.0000376718e"2^(13.117 - 8.74678y + y'^f 
W^^4y) = e"22'(-0.454214 + 0.759976y - 0.432314y2 + 0.0971138y^ - 0.00730079/ 

+ 0.0000251146/) 

(C.14) 

Iron (5®Fe) 

W^{y) = e"22'(62.3888 - 160.428y + 152.644/ - 67.2779/ + 14.478/ - 1.43665/ 

+ 0.0525291/) 

l+/(y) = e“22'(0.318309 - 1.27323y + 1.99188/ - 1.54562/ + 0.622264/ - 0.122277/ 

+ 0.00921525/) 

W}^{y) = e“2*'(-4.45633 + 14.6422y - 18.2579/ + 10.8919/ - 3.2296/ + 0.446836/ 

- 0.0220016/) 

W^iy) = e"22'(-4.45633 + 14.6422y - 18.2579/ + 10.8919/ - 3.2296/ + 0.446836/ 

- 0.0220016/) 

1+00 (y) ^ £-22/(4.22872 - 6.76595y + 3.79067/ - 0.867433/ + 0.069506/) 
l+/,(y) = e“22'(0.143378 - 0.229404y + 0.144606/ - 0.0422756/ + 0.00486921/) 
l+/,(y) = e“22'(-0.778655 + 1.24585?/ - 0.741661/ + 0.194658/ - 0.0183967/) 

1+01 (y) ^ 1.24585?/ - 0.741661/ + 0.194658/ - 0.0183967/) 

l+^^„(y) = e“22'(-16.2427 + 33.8776?/ - 25.2342/ + 8.30471/ - 1.20334/ + 0.0604243/) 
l+/^„(y) = e"22'(-0.213631 + 0.598168?/ - 0.622338/ + 0.308014/ - 0.0735211/ 

+ 0.00669858?/^) 

WW^„(y) = e"22'(l.16019 - 3.24853?/ + 3.31473/ - 1.54264/ + 0.325833/ - 0.0253084/) 
W^^4y) = e“22'(2.99085 - 6.23805?/ + 4.81422/ - 1.74483/ + 0.288128/ - 0.015993/) 

(C.15) 

Nickel (58Ni) 

l+“(y) = e"2^(66.9246 - 175.389?/ + 169.877/ - 76.127/ + 16.6597/ - 1.6839/ 

+ 0.0628067/) 

l+/(y) = (0.0795762 - 0.318305?/ + 0.548985/ - 0.503018/ + 0.250492/ 

- 0.0603789/ + 0.00545169/) 

l+]0(y) = e“2^(-2.30773 + 7.63937?/ - 10.3404/ + 6.95311/ - 2.30652/ + 0.350525?/^ 

- 0.0185041/) 

W^{y) = e"2^(-2.30773 + 7.63937?/ - 10.3404/ + 6.95311/ - 2.30652/ + 0.350525/ 

- 0.0185041/) 
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= e"2^(5.4697 - 8.75152y + 4.8845V - 1.10715y^ + 0.0875404/) 

Wll{y) = (0.00977975 - 0.0156476y + 0.0136707/ - 0.00592935/ + 0.00140426/) 

Wl^,{y) = e"2^(-0.231284 + 0.370054y - 0.264922/ + 0.0935201/ - 0.0110873/) 
Wll{y) = e“2^(-0.231284 + 0.370054y - 0.264922/ + 0.0935201/ - 0.0110873/) 
W^^„{y) = e“2^(-19.1326 + 40.3764?/ - 30.3339/ + 10.0435/ - 1.4629/ + 0.0741493/) 
W}}^„{y) = e“2^(-0.0278969 + 0.0781112?/ - 0.0956406/ + 0.0607914/ - 0.0211633/ 

+ 0.00276687?/^) 

W}^^„{y) = e"2^(0.659741 - 1.84727y + 2.0953/ - 1.10461/ + 0.25912/ - 0.0218459/) 
W^^„{y) = e"2^(0.809015 - 1.7073?/ + 1.48687/ - 0.692274/ + 0.145722/ 

-0.00939131/) (C.16) 
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